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Abstract 

Electric and magnetic fields created by moving uniform lattice of Josephson vortices 
in the magnetic field parallel to the layers are calculated in the frame of exactly solvable 
model. At large velocities of the vortex lattice the plasma oscillations of superconducting 
electrons are excited by vortex motion, this results in interference features in I-V curves at 
low temperatures. The spectrum of electromagnetic radiation by moving vortices contains 
peaks related to the excitation of plasmons and to the Cherenkov radiation. 

It is known that in high-Tc superconductors there exist low-damping plasma oscillations f^- 
[0 in the range of hundreds megahertz - terahertz. The presence of an underdamped eigenmode 
leads to a resonance behavior of the driven modes when frequency and wave- vector of the driven 
force are close to those of the eigenmode. In particular, frequencies of the electromagnetic field 
of sliding lattice of Josephson vortices created by magnetic field parallel to the layers can 
enter the range of plasma modes if the lattice is moving fast enough. The velocity of such 
vortices driven by transport current perpendicular to the layers may be very large, because the 
amplitude of the order parameter A inside Josephson vortices is perturbed very slightly. Thus, 
at large enough voltages across the superconductor the resonance conditions can be satisfied. 
This results in the generation of plasma oscillations and influences the I-V curves. Below we 
study the motion of uniform lattice of Josephson vortices in an infinite crystal in magnetic field 
which is small enough, so that non-linear cores of the vortices do not overlap. We assume that 
the characteristic frequencies of the problem are small in comparison to A. 

In the presence of vortices the spectrum of the eigenmodes is different from the zero-field 
spectrum due to oscillations of vortex displacements leading to the sound-like |IT|,[0)|@| spectrum 
of the modes. In contrast, the response of a superconductor to the motion of the vortex lattice 
driven by the transport current calculated in our study, is not related to the oscillations of the 
vortices, driven modes being related to the plasmons with the plasma edge. 
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The current distribution in the vortices is influenced by their motion, this makes calcula- 
tions of the vortex dynamics difficult since usual perturbational approach is not valid at large 
velocities of the lattice. To solve the problem strictly one must find solution of Maxwell's 
equations with expressions for charge and current densities in layered superconductors. Having 
in mind high-Tc superconductors we shall consider the order parameter having the symmetry 
close to d-pairing. Expressions for current and charge densities are non-linear functions of the 
phase differences between adjacent layers, which are rather complicated functions in general 
case. Quasi-particle current densities depend on phase differences, on the supecronducting mo- 
mentum Pn = {l/2)iqxn — (l/c)A„ and, in addition, on the gauge-invariant scalar potential 
fin = (l/2)(9xn/<9t) + $n, where A„ is vector potential, is electric potential, and Xn is 
the phase of the order parameter in layer n. (We assume here h = l,e = 1). In general case 
potential fin to be found from the Poisson's equation, however in the vortex problem at low 
temperatures considered here, T <C A, the quasi-particle branch imbalance, and, hence, poten- 
tial fin can be neglected provided (ro/d)'^ <C 1, where tq is Thomas- Fermi screening radius, and 
d is the period of the crystal in the direction perpendicular to the layers. Current density along 
the layers can be described by expression for linear response since the characteristic value of the 
current is determined by the critical current density jc in the perpendicular direction, which is 
small in comparison to the critical current in the parallel direction. Then the expression for 
current density in layer n, presented by means of Fourier transformation has the form 
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The first term here describes the superconducting current, and the second is related to the 
quasi-particle current, —iu and zq correspond to the time derivative and to the gradient in the 
direction parallel to the layers. The frequency dependence of conductivity cry (a;) is determined 
by symmetry of the order parameter and depends on the momentum scattering time. We shall 
use the expressions from Ref. to describe cr||(ti;). 



The density of the superconducting current between layers n and n + 1 is given by 

where cpn is the gauge invariant phase difference. We shall adopt an exactly solvable model 
substituting sin by the saw-tooth function 

Jxlifn) = jc arcsin sin (pn, (1) 
as it was done in [0. Such a model was used to study the Josephson vortex flow in the 



linear approximation on velocity in |]r2|. Following (|I|) we shall use linear relation between the 
quasi-particle current and the phase difference 

jff = -^uJa^{co)^n/{2d). 
The expressions for the current densities we insert to the Maxwell equation 
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in which we shall use discrete approximation for spatial derivatives in z direction perpendicular 
to the layers, and express magnetic field in y direction as and P^i 
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where x is the coordinate in the direction of the vortex fiow. 

In the displacement current we take into account only the component along z axis, since 
due to the strong anisotropy of a layered superconductor the plasma frequency in the direction 
parallel to the layers, flp = c/X, is much larger, than typical frequencies of the problem, 
which are of the order of the plasma frequency for perpendicular direction, Up = Cz/Xc- Here 
Cz = c/ -y/e, e is a dielectric constant in z direction, A and Ac = c/ ^/87!^djc are the screening 
lengths for the current flowing parallel and perpendicular to the layers, respectively. Then we 
get the equations 



The solution of these equations yields spatial and temporal dependencies of electric and mag- 
netic fields. 

Using the model current-phase relation ([|) we can find the exact solution of equations 
by means of the Fourier transformation. In principle, we must take into account a bending of 
the vortices due to their motion. However, we shall limit our study here by the case when the 
displacement of the vortex center xo{y) along x axis is small at distances y ~ Ac. Estimation 
of the vortex bending based on equations of a balance between forces acting to the vortices 
and created by the transport current and by the vortex currents shows that the deformation of 
vortices can be neglected in magnetic fields H ^ Hd provided 

Ac TTd\n{X/d) jcH 

We assume that H is large in comparison to Hd and condition is satisfied. 
Then for the triangle lattice sliding with velocity u we get 

{ul -uj^ + iujujr)f - 2c^qKP/t = U6{q - uo/u) (6) 
-{c^qK/2d^)ip + ~uj^ + iujnr)P = 0. (7) 

Here K = 2sin/i;/2, < vr is the wave number obtained from the discrete Fourier transforma- 
tion with respect to the layer number, uj^ = Up{l + A^g^), Qp = Qp{l + X'^K'^/d'^), Ur = 4:'n-a±e, 
Qr = 4:7ca\\ and, finally. 



H = A7T^{ujl/iuj) exp + m/2)qX - imkZ/d]. 

l,m 
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Solution of equations (^0) for the phase differences takes the form 

Lp = Ii{l + X^K^/£)5{q-u/u)/D, (8) 

where D is the determinant of equations (1^-0). Zeros of D yield the spectrum of the eigenmodes, 
i. e. of the plasmons. At high frequencies and low temperatures, when frequencies of dielectric 
relaxation Reuj^ ReQr, are small enough, we get 

D = [{ul - uj^){l + X^K^d^) + ulXlq^ - iuj[ujr{l + X^K^/d^) + ^.^^ - uo^ + uolXlq^)/^^]. 

Variation of the solution with frequency increasing can be illustrated by the example of the 
component of the electric field = —i{uj/2d)(f. At frequency u = qu {q = 27t/X) the slowly 
varying along z part of the field Ez created by a single vortex has the form 
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(9) 

The total field is described by the sum of the fields created by individual vortices. One can see 
that at small frequencies the decaying length is of the order of A. When frequency becomes 
larger, than the plasma frequency, the real part of the exponent in becomes small. This 
means that in the range of plasma oscillations the decaying length for the electric field is large 
and is determined by the damping of plasmons. 

In order to find the relation between transport current jtr and the velocity of the vortex 
lattice M, we multiply equation (j^) by {c^ /X^^)d(pnldx, and @ by idPn/dx, then we add these 
equation, integrate them by x and sum up over all n. Then in the left-hand part of the resulting 
equation we get an expression which is a full derivative over x, the integral of which can be 
reduced to —2d'^ J2n Hy This expression is related to the force acting on the vortex lattice 
and is proportional to the half sum of the magnetic fields at the opposite sides of the sample 
{i. e. to the externally applied magnetic field), and to the difference of the magnetic fields, 
which by means of equation (H) can be expressed in terms of the transport current. The right- 
hand side we express by means of the Fourier transformation using solutions of equations (^0). 
Finally, we get 

IQT^^djtr = J dqdkuulleuril + X^K^/d^f + ^IXK\^ / d^)]/\D\''Y.e'^^+'^/^^^+'"'^^''^ . (10) 

l,m 

Integrating (p!0|) we find a relation of the transport current to the velocity of the vortex lattice, 
and, hence, to the average electric field E = 2imLU (where rti is the vortex density). 

In the limit of small velocities, u ^ {d/X)c, one can neglect the effect of velocity u on the 
shape of vortices. In this limit I-V curve is described by the Ohm's law with the resistivity 

1 2nLdXj 
a_L + cr||(A/Ac)2 n 

where Xj = dXd X. The last factor in (|l^) describes the relative volume of the superconductor 
occupied by the non-linear region. One can see an analogy to the Bardeen-Stephen law with 



4 








0.02 0.04 0.06 0.08 0.1 0.12 0.14 



Figure 1: Typical relation between normalized transport current and vortex velocity in mag 
netic field H = 20H^i. 



the area of the Abrikosov vortex core substituted by the area of the non-linear part of the 
Josephson vortex dXj. 

An analytic expression for the dependence of u on jtr may be obtained easily also for 
frequencies of the order of the plasma frequency for the case of small plasmon damping 



Jtr = 3c 



sinh h 



cosh h — cos a 



IjJr OJ. 
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The last factor in this equation describes the oscillations at the I-V curve appearing as a result 
of interference of the plasmons created by different vortices. In the limit of large velocities the 
I-V curve approaches the Ohm's law again, but with resistivity p = 2tc ^ nidXj / a{ujp) . The 
shape of the I-V curves is sensitive to the magnitude and to the anisotropy of the quasi-particle 
conductivity. An increase of the conductivity results in an increase of the plasmon damping 
and in a decrease of the peaks. The regime of very small damping is easily achieved in the 
case of isotropic pairing when the order parameter has no nodes and the quasi-particle density 
is exponentially small. In the case of d-pairing the damping is higher and peaks at the I-V 
curves are smaller. In Fig.l we present the dependence jtr{u) calculated numerically under 
assumptions that the real part of the conductivity is described by expressions from [|lOl and 
is decreasing at frequencies above the quasi-particle momentum scattering rate as Note 
that the condition (H) is difficult to satisfy at large currents, therefore the shape of the I-V 
curve near the maxima in Fig.l may need more detailed calculations which take into account 
the bending of the vortices due to their motion. 

Thus, in the range of plasma oscillations the regions of negative differential conductivity 
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Figure 2: Typical dependence of energy flow radiated at the frequency of the first harmonic, 
H = 20H^i. 



appear in which the uniform flow of the vortex lattice is unstable. These regions are separated 
by the stability regions. At large voltages the uniform motion is stable. 

Excitation of plasma oscillations is related to an emission of electromagnetic radiation. We 
calculate the energy flow traveling along x axis at frequency un = 2ttuN/X, where is the 
harmonic's number. To do this we find the Pointing vector using expressions for electric and 
magnetic fields calculated by means of (Bf). 
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It is evident that the energy flow increases at resonance frequencies, when the denominator 
ID I is small. In the case of triangle lattice we find that the energy of odd harmonics contains 
maxima in the region of plasma frequency, while even harmonics has additional sharp peaks at 
the velocities of the lattice close to the velocity of the light in the superconductor Cz playing 
a role of the Swihart velocity of the tunnel junction. The latter peaks correspond to the 
Cherenkov radiation. The energy flow created by the first harmonic, S{u), is shown in Fig. 2, 
where the energy flow is presented in units of 5*0 = h^(?ujp/ [l&e^X^Xc). Using values typical 
for BSCCO, e = 25, A = 0,2 /im, = 60 /xm, we find 5*0 ~ 10 W/cm^. The shape of the 
curve and the peak's size strongly depend on the conductivity value and on the anisotropy of 
the conductivity. The values of the lattice velocity corresponding to the maximums at S{u) do 
not necessarily correlate with the peaks of the conductivity, thus, the peaks of the radiation in 
the range of the plasma frequency can be observed in the regions of the stability of the uniform 
vortex lattice flow. 
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The peak in the energy flow at u > Cz is much higher than those near the plasma frequency, 
it can be approximated by the expression 

SJl 

47r2{[l + /iAr2(l - mVc2)]2 + hN^u^iu = c,)}' 

where h = An'^XXcnL. 

Note that the frequencies of the radiation at the vortex motion with velocities close to 
belong to the frequency range below the superconducting gap value only if the magnetic field 
is small enough, i. e. when the period of the vortex lattice X is comparable to Ac- In larger 
fields the corresponding frequency approaches A. However, this must not affect the results 
considerably, because at such frequencies the part of the response corresponding to damping is 
very small in pure, tT^ ^ 1, superconductors. 

Note also that the effects considered here are related to the range of higher voltages and 
frequencies, than those typically used in experimental studies of Josephson vortices, e. g. in 
[0, in which a non- Josephson emission from the moving vortices was reported in BSCCO. 
As long as we consider here the uniform motion of the lattice in the infinite crystal, such a 
radiation in our treatment is absent. 
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